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We discuss the effects of disorder in time-reversal invariant topological insulators and supercon- 
ductors in three spatial dimensions. For three-dimensional topological insulator in symplectic (All) 
symmetry class, the phase diagram in the presence of disorder and a mass term, which drives a 
transition between trivial and topological insulator phases, is computed numerically by the transfer 
matrix method. The numerics is supplemented by a field theory analysis (the large-A^/ expansion 
where Nf is the number of valleys or Dirac cones), from which we obtain the correlation length 
exponent, and several anomalous dimensions at a non-trivial critical point separating a metallic 
phase and a Dirac semi-metal. A similar field theory approach is developed for disorder-driven 
transitions in symmetry class AIII, CI, and DHL For these three symmetry classes, where topologi- 
cal superconductors are characterized by integer topological invariant, a complementary description 
is given in terms of the non-linear sigma model supplemented with a topological term which is a 
three-dimensional analogue of the Pruisken term in the integer quantum Hall effect. 



I. INTRODUCTION 

A topological phase is a gapped quantum state of mat- 
ter with non-trivial quantum correlation (entanglement). 
Their key distinction from trivial states a band insula- 
tor, say, lies in the fact that they support stable gapless 
boundary modes when terminated by a boundary. A sys- 
tematic classification of topological insulators and supeiu 
conductors in all Altland-Zirnbauer symmetry classesEJ 
and in all dimensions was obtained, revealing mod 2 and 8. 
periodicities in dimensionality and symmetry classes.Bcl 
It was found in any dimensions there exist topologically 
non-trivial insulators or superconductors for five out of 
ten Altland and Zirnbauer symmetry classes. 

In three spatial dimensions, among the five symme- 
try classes is symplectic (spin-orbit) symmetry class 
(class AJJXri in which Z2 topological insulators can be 
realized.Lrii3 A number of materials with non-trivial Z2 
topological number have been discovered, starting from 
Bismuth- Antimony alloys. (See, for a partial list of such 
materials, Refs. ^l|-jl^) 

Topological superconductors that are possible in three 
dimensions are all time-reversal symmetric; they are real- 
ized in symmetry class CI (conserved SU(2) spin rotation 
symmetry), AIII (partially conserved SU(2) spin rotation 
symmetry), and Dili (no SU(2) spin rotation symmetry). 
In these symmetry classes topologically different gapped 
states are characterized in the bulk by an integral, as 
opposed to Z2, topological number u defined in terms 
of the spectral projector in momentum space. In these 
topological superconductors, the Bogoliubov quasiparti- 
cles are fully gapped in the bulk by the (mean field) pair- 
ing gap, and the system is a thermal insulator. (When 
z-component of spin is a good quantum number as in 
symmetry class CI and AIII, spin transport can be dis- 
cussed; topological superconductors in these classes are 



a spin insulator in the bulk.) On the other hand, when 
such three-dimensional (3d) topological superconductors 
are terminated by a 2d surface, they can support gapless 
and stable Dirac (class AIII and CI) or Majorana (class 
Dili) surface modes, the multiplicity (number of Dirac or 
Majorana cones) of which is specified by the bulk topo- 
logical number u. The fermionic quasiparticles in the B 
phase of '^He (the BW state) realize an example of topo- 
logical superconductor (superfluid) in class DHL Singlet 
BCS pairing models on the diamond and cubic lattices 
that realize a topological superconductor in class CI have 
also been proposedJlj'ES 

While topological insulators and superconductors are 
stable against weak perturbations, with strong perturba- 
tions, they can undergo a (continuous) quantum phase 
transition into a different phase. For example, if one can 
tune the chemical potential in the BW state, there is a 
transition between a trivial phase (strong pairing phase) 
and a topological phase (weak pairing phase). Likewise 
in the quantum Hall effect (QHE), strong enough disor- 
der can also destroy the topological insulators and su- 
perconductors, and drive a system, through a continuous 
quantum phase transition, into a metal or other topolog- 
ical/ trivial phases. 

This paper is devoted to study such quantum phase 
transitions among topological (electrical, thermal or 
spin) insulators, trivial insulators, and a metal, in sym- 
metry class All, and in three Bogoliubov-de Gennes 
(BdG) symmetry classes CI, AIII and DHL These quan- 
tum phase transitions can be induced by changing param- 
eters that modify the band structure and wavefunctions, 
and also by introducing disorder; both of these perturba- 
tions are assumed to be preserving discrete symmetries 
defining these symmetry classes. 
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A. main results and outline of the paper 



II. DISORDER IN THREE-DIMENSIONAL 
TOPOLOGICAL INSULATOR IN SYMMETRY 
CLASS AH 



a. symplectic (All) symmetry class Wc start, in 
Sec. |l[ by introducing a 3d Dirac-type Hamiltonian on 
the cubic lattice in symplectic (All) symmetry class, 
which realizes a topological insulator, a trivial insulator, 
and a phase transition between these two. In Subsec. 
II B, the phase diagram in the presence of a mass term 
and disorder is determined numerically by the transfer 
matrix method. (For earlier studies, sec, for example, 
Refs. 0-|o[) Wc observe that with disorder the phase 
boundary between topological and trivial phases is renor- 
malized, and in particular the topological phase gets "en- 
larged" . This means, as we pass through along a particu- 
lar cut in the phase diagram, we can turn a trivial insula- 
tor into a topological one by increasing disorder strength 
("topological Anderson insulator"). 



The numerical analysis is supplemented with a field 
theory analysis (the large-iV/ expansion of a Gross- 
Neveu-type model where Nf is the number of valleys or 
Dirac cones), in Subsec. II C, from which we obtain the 
correlation length exponent, several anomalous dimen- 
sions at non-trivial critical point, and dynamical expo- 
nent. 



b. BdG symmetry classes AIII, Dili and CI In 
time-reversal symmetric superconducting symmetry 
classes (class AIII, Dili and CI), we will develop, in Sec. 



Ill, a similar large- iV/ field theory approach in terms of 
chiral Gross-Nevcu- or Nambu-Jona-Lasinio-type models 
and their large- TV/ expansion. 

While description in terms of these fermionic field the- 
ories are reliable for a region of the phase diagram where 
spin or thermal conductivity is relatively small (in par- 
ticular near the clean fixed point separating trivial and 
topological superconductors), deep in a spin or thermal 
metallic phase, on the other hand, a complementary de- 
scription in terms of t he no n-linear sigma model (NLcrM) 
is possible [Subsec. |IIID |. For these three symmetry 
classes, where topological superconductors are character- 
ized by the integer topological invariant v, the NLctMs 
can be supplemented with a topological term which is 
of integer type (related to ^^{0 / H) = Z of the NLctM 
target space G/H); it is a higher diiiacjisional analogue 
of the Pruiskcn term in the QHEE2rc3. We will com- 
pute the theta angle in the topological term in terms of 
microscopic parameters (the fermion mass and the imag- 
inary part of the fermion self-energy in the self-consistent 
Born approximation). While we do not have an analyti- 
cal tool to analyze the effects of the topological term on 
the nature of Anderson metal-insulator transitions, we 
will speculate on possible RG fiows in terms of diagonal 
spin or thermal conductivity and the theta parameter. 



A. Hamiltonian and symmetries 

stapti 
lianESt 



Let us staf±,|by considering the following tight-binding 
Hamiltoniai 



X fe=l,2,3 



t ^ t ^ , 



(w + 3r)E4/3cx, 



(1) 



where c]^, c^ represents the four-component fermion cre- 
ation/annihilation operator defined on a site x on the 
3d cubic lattice spanned by three orthogonal unit vec- 
tors ek=x,y,z, and the 4x4 gamma matrices in the Dirac 
representation are given by 



cjk 
Uk 



,/3 



1 

-1 



75 



1 

1 



, (2) 



with i7fc=i_2,3 being the standard Pauli matrices. In 
momentum space, the Hamiltonian is written as Hq = 
J2k4j^oO^)Ci^, where 'Ho(k) is given by 

3 

Ho(k) ==E«'«(^K + ^4(k)/3, 

a=l 



with 



daCk) = t sin ka, (a = 1,2, 3), 
3 

(i4(k) = m + ^^(1 — cos ka), 



(3) 



where k = {k^,ky,kz) G [—tt^tt)^. For arbitrary values 
of the parameters t,m,r £ M, the Hamiltonian is time- 
reversal symmetric, («tTy)'H5(— k)(— icTj,) = Ho(k). In the 
following we fix the Wilson parameter (r) and the hop- 
ping amplitude (t) as r = 1 and t = 2. 

As we change m in the Hamiltonian, we can realize the 
trivial insulator, the Z2 topological insulator, and also a 
quantum critical point separating trivial and topological 
insulators: 

• For m < —6, —4 < m < —2 and < m, i?o realizes 
a trivial insulator whereas 

• for —6 < m < —4 and —2 < to < 0, i^o realizes a 
topological insulator. 

In the following we will mainly focus on the region near 
TO = 0. 

We add, to Hq, the following disorder potential 



(4) 



which is time-reversal symmetric, and where the random 
variable Ux is distributed between —W/2 and W/2 ac- 
cording to the box-distribution. 
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When we are close to a quantum critical point, say, at 
m = 0, the tight binding Hamiltonian can be described 
in terms of its continuum limit, which takes the form of 
3d massive Dirac Hamiltonian 



H = ■Ho + V, Ho = -idkctk + mP, 



(5) 



where the summation over repeated indices k — 1, 2,3 is 
implicit, and V represents a disorder potential. For any 
realization of disorder, the continuum Hamiltonian (|^) 
satisfies time-reversal symmetry 



i<TyH* {—i<Ty) = Ti. 



(6) 



The clean part of the Hamiltonian Ho realizes, de- 
pending on the sign of the mass term m G M, a topolog- 
ically trivial and non-trivial insulator in symmetry class 
All, each of which is characterized by the vanishing and 
non- vanishing, respectively, of the Z2 invariant. When 
m = and in the absence of disorder, the spectrum is 
gapless whilst the dc conductivity is zero, (Jxx = 0. This 
should be compared with the "uniiiersal" finite conduc- 
tivity in the 2d Dirac Hamiltonian£3 



B. 



numerics 



c. transfer matrix method In this section, we will 
study numerically effects of disorder in the 3d Dirac 
model. Below, we use the transfer matrix to study the 
random Hamiltonian H = Hq + V. To this end, we write 
the "Dirac equation" for the single-particle wavefunction 
"(/'(x, ct) = ip(r,a)n at a given energy E in the following 
form: 



1 -t2 
1 



l-t2 



{(3 + ita,)[n^"^ - E]i:^ 
{P + ia^){l3 -iaz)i^n-i- (7) 



Here we decomposed the 3d spatial coordinates x in 
terms of its x- and y- components r and its z-component 
n, X = (r,n), and "H^"-* is the Hamiltonian for the rtth 
cross section. As shown in Appendix ^ this can be rep- 
resented as 



1 oy 



where = + - -^1 



(8) 



The size of each slice with fixed n is i x i, so that A/^") 
is a 4i X AL matrix. We will denote the length of the 
z-direction by N . 

We study numerically the dependence of the small- 
est Lyapunov exponent of the transfer matrix M, as a 
function of the width of the quasi-one dimensional 
system. Lyapunov exponents are self-averaging random 
variables for an infinitely long quasi-one dimensional sys- 
tem, N ^ 00. 



The eigenvalues of the Hermitian matrix M^M are 
written as exp{zL2Xj) with < ATi < X2 < . . . < 
The decay length is then given by 



N- 



NXi 



(9) 



The decay length is a finite and self-averaging length 
scale that controls the exponential decay of the Landauer 
conductance for any fixed width L of the infinitely long 
quasi-one dimensional system. It is of course impossi- 
ble to study infinitely long quasi-one dimensional lattice 
models numerically and we shall approximate with 
^ obtained from the Lyapunov exponents of a finite 
but long quasi-one dimensional network model made of N 
slices. In our numerics we have set iV = 1^ 10^ ^ 5 x 10^. 

As shown by MacKinnon and Kramer ,lZI criticality can 
be accessed from the dependence of the normalized decay 
length 



(10) 



on the width L of the quasi-one dimensional system. As 
usual, Al is required to be written as A^ = f{L/^) on the 
basis of the single-parameter scaling hypothesis, where ^ 
is the (3d) localization length. At the critical point the 
localization length diverges according to the power low 
^ ^ |a; — Xc\~^^'^ upon tuning of a single microscopic 
parameter x (disorder strength W or mass m in the fol- 
lowing) . 

The numerically computed normalized decay length is 
shown for to = and for different disorder strength in 
Fig. I]. The plots for different values of L (=6, 8, 10, 
and 12) cross at a single point W ^ 3.9, signaling an 
insulator-to-metal transition. The inset shows Al as a 
function of L/^ with the help of the one-parameter scal- 
ing ansatz. In a metallic region fits well with a typical 
metal-insulator transition behavior, while in the insulat- 
ing phase L-dependence gets severer at weak disorder. 
This is so since in the insulating phase, the presence of 
the band gap would complicate the scaling analysis. The 
fact that in the insulating side the data are not well fit by 
the scaling curve also makes the unambiguous determina- 
tion of the critical exponent difficult. A similar problem 
was observed in the metal-iuaulator transition in the 2d 
quantum spin Hall systems .Ej^3 A special effort on our 
3d model will be made elsewhere. 

d. structure of the phase diagram In Fig. ^, the nu- 
merically determined phase diagram is shown. Observe 
that as we increase disorder the phase boundary separat- 
ing trivial and topological insulators gets shifted toward 
the large mass region, i.e., the trivial insulator region 
becomes topologically non-trivial as disorder strength 
increases. This disorder-induced topological iasulating 
phase has been found first ip-lwo dimensionsE2rc3 and 
recently in three dimensionsE3Ej, and thereby referred 
to as "topological Anderson insulator". (However, one 
should note that the determination of the phase bound- 
ary between the two insulating phases is rather difficult 
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FIG. 1. Normalized decay length Al = A/L as a function 
of disorder strength W for m = 0.0. Inset: A fit of the data 
shown in the main panel with the help of the one-parameter 
scaling ansatz in the vicinity of the critical point. 
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FIG. 2. The numerically obtained phase diagram in the 
{W, m) plane. Lines are guide for eyes. 



numerically since in an insulator the L dependence of the 
decay length is intrinsically small.) While it is not de- 
picted here, in the strong disorder limit, we expect there 
should be another phase boundary separating the metal- 
lic phase and an insulator. 



C. field theory 

Wc now study the metal insulator transition from field 
theory point of view. Our starting point is the continuum 
limit of the lattice tight-binding model, Eq. (0), 



(11) 



With time reversal symmetry, the disorder potential V 
can be expanded in terms of ax,zTy, ctoTo,x,z, and ayTy, 
and hence parameterized by six real parameters charac- 
terizing disorder. Below, we assume, for simplicity, dis- 



order potentials are distributed according to a Gaussian 
distribution, and their variances are all equal. 



P{V) oc exp 



--Ltr4VVt 
25 



(12) 



where 5 is a parameter that characterizes the strength of 
disorder. 



1. fermionic replica method 

Wc will use the fermionic replica method for quenched 
disorder averaging. In the fermionic replica method, the 
discrete symmetry of the problem, time-reversal symme- 
try, is implemented as a continuous symmetry: invariance 
under 0{2Nr) x 0(2AV) rotations in the replica space 
where is the number of replicas. 

In order to study Anderson localization physics, we 
are after the properties of the single particle retarded 
(R) and advanced (A) Green's function, (±iO+ — H)^^, 
and in particular, the product thereof, where 0^ > 
is a small imaginary part of the energy needed to regu- 
larize poles. For a product of the retarded (R) and ad- 
vanced (A) Greens functions, cx (iO+ — H)^^(iO+ + H)^^, 
its generating function can be expressed as a fermionic 
functional integral, Z ~ J cxp(— J cfixC), where 



£ = * (m + O+A) 



(13) 



Here, the total number of components of the fermion field 
^, * is 



(Dirac) (g) (TR) ® (R/A) (g) (replica) 
^ (T(g}T(g} (R/A) (E) (replica) = 16Nr 



(14) 



where A is a diagonal matrix which is ±1 for the re- 
tarded/advanced sector, respectively. The fermionic field 
^' incorporates the time-reversal symmetry grading, and 
can be written as 



* = 



1 



.*^4( 



X 



V2 V ^^<^yX 



(15) 



where x, x are two-independent fermionic fields (we have 
suppressed the replica indices.); the subscript (• • • )t- 
means the block structure displayed in Eq. (|l^) repre- 
sents the time-reversal symmetry grading. It will prove 
convenient to rewrite Eq. (^5|) in the spin grading in 
which the block structure displays subspaces with cr = ±: 



where 



I zy 



K := (xt^Xi), 7 := (Xi,-Xt)' 

{jz ~ Ty)/^. 



(16) 



(17) 
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With these, the action can be written correspondingly as 
* {in + 77 A) * 



By introducing four component spinors by 



, = i>^C'^ = ^/'lay, (19) 



we end up with the functional integral Z — 
J 'Dlijj, ip] exp (— J (fix £) with the Lagrangian 



(20) 



,.fa=l 



When 0+ — 0, the action is invariant under 0(ANr) ro- 
tations, K — >■ Uk, 7 — > J77, U G 0(4A'^r)- In the metallic 
phase where a small imaginary part in the self-energy is 
generated spontaneously, this symmetry is broken down 
to 0{2Nr) X 0{2Nr). The resuhing Nambu-Goldstone 
mode is the diffuson in the metallic phase. 

We can now perform the disorder averaging accord- 
ing to P{V), resulting in the generating functional Z = 
J 'D[ip, ip] exp (— J cPx C) with the Lagrangian 



„b=l 



(21) 



Whcn|-0."'" is zero, this is a 3d version of the Gross-Neveu 
modeO with 0(4iVr) internal and 0{Nf) flavor symme- 
tries. 



2. saddle point 

The disorder-induced four-fermion "interaction" can 
be decoupled by an auxiliary matrix field Wab leading 
to Z = / V[^P, iP] J V[W] exp (- / (Px C) , 

Nf iNr 

/: = ^ 51 i^'^O^c.b + 0+A,6 + Wab) ^,b 
L=l a,b=l 

+ [^W^] ■ (22) 

Here, we have introduced TV/ flavors of fermion fields 
■010 (t = 1,---,A^/)- Equation (^ ) reduces to the sit- 
uation we have been discussing when Nf = 1, whereas 
Nf > 1 corresponds to the case with Nf "valleys" . The 
replicated Lagrangian ( ]2^ ) can be derived from a model 
similar to ( pT| ) once we assume all the flavors are coupled 
equally by a disorder potential. (See Appendix ^) We 
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FIG. 3. The mean field (self-consistent Born) phase diagram 
near the Dirac point. Arrows indicate the renormalization 
group flow. 



will discuss the four-fermion "interaction" term in terms 
of the large Nf expansion below. 

As a first step, we now look for a spatially homoge- 
neous, replica symmetric saddle point solution, by set- 
ting 



Wab = vKb- 



(23) 



This reduces the symmetry 0{4,Nr) 0{2Nr) x 0{2Nr). 
The self-consistency conditions for 77 is given by 



1 



(27r)3 fc2 + m2 + r;2 ANfg 



0, 



(24) 



where we have introduced the ultra-violet cutoff Auv 
Anticipating spontaneous breaking of 0{ANr) symmetry, 
we have set 0+ = 0, as its sole role is to select the pattern 
of the residual symmetry, 0{2Nj.) x 0(2A'^r)- 

The resulting mean-field (or self-consistent Born) 
phase diagram is given in Fig. ^. We identify the phase 
with ?y 7^ (the ordered phase with spontaneous 0{ANr) 
symmetry breaking) as a (diffusive) metallic phase. This 
metallic phase appears because disorder creates states 
that fill the pseudo band gap of the Dirac fermions. On 
the other hand, there are two, trivial and topological, 
insulating phases with 77 = which are separated by a 
phase boundary which runs vertically at 771 = 0. 

In the lattice Dirac model (^, a unitary transforma- 



tion Cx — «75Cx 



-c^i75 relates two Hamiltoni- 



ans, one with the parameter (m, r) and the other with 
(— 777,, — r); i.e., Hamiltonians with the mass jti and —to 
arc not unitary equivalent if r is held fixed. On the other 
hand, the continuum Dirac model ([ll| ) and its replicated 
counterpart (]2^) are invariant under a Z2 transformation 
"00 — >■ 7500, i^a — >■ "0075, 777 — >■ —TO. For this reason, while 
the numerical phase diagram Fig. (|2|) is asymmetric un- 
der TO -!->■ —TO, the mean field phase diagram Fig. ^ is 
symmetric under m -H- —to; the "topological Anderson 
insulator" phase is not possible in the continuum model 
with linear dispersion. 
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3. 1 / N f -expansion 

Just at the Dirac point, the disorder, which appears as 
a four fermion interaction in the generating functional, 
is irrelevant, whereas the mass term is relevant, from 
the power-counting. The transition between trivial and 
topological insulators is then described by the clean Dirac 
point. We will now turn our attention to the non-trivial 
fixed point at (rn,g) — (0, .gc), where 7^ is the non- 
trivial solution to the saddle point equation Eq. . To 
include 1 /Nf corrections systematically, we consider the 
Lagrangian 

Z^Z^^^,aWati^,b + \ZwWabWab. (25) 



where repeated flavor and replica indices are implicitly 
summed, and we have introduced Euclidean gamma ma- 
trices 7fc = Qffe , g is the bare coupling and Z^ , Zw are 
wavefunction renormalization constants; we have also 
rescaled the bosonic field Wab properly. These renormal- 
ization constants must be adjusted at each order in l/iVy 
to keep vertex functions (in particular, the fermion prop- 
agator, the fermion-boson vertex, and the boson propa- 



gator) finite. (See, for example, Refs. 36-39 



To leading order in l/Nj, as we have seen in the sad- 
dle point analysis, the expectation value of W (=77) is 
spontaneously generated ii g > gc- the self-consistency 
equation is 



(fik 



1 



(27r)3 + r/2 NfD^g 



(26) 



where D^{= 4) is the dimensionality of the Euclidean 
gamma matrices. The expectation value 77 is a physical 
observable which sets a scale, and hence independent of 
the UV cutoff. We will go to the non-trivial branch g > gc 
of the saddle point equation ( p6[ ) for the most part of our 
analysis below. This equation relates the bare coupling 
constant g to the UV cutoff and a physical scale t], and 
can be considered as a renormalization condition. 
The fermion propagator to leading order is given by 

(V'm(fc)'0K6(fc)) = S.i^Go.abik), 

g^\k) = 2Z^{ika^ + ^^ab), (27) 

where we can choose Z^ = 1 and E = 77 to this order. The 
factor of 2 reflects the Majorana nature of the fermion 
fields. 

The boson two point function to leading order can be 
computed by summing over fermion bubble diagrams (see 
Fig. I), 



{Wab{q)Wcd{-q')) 



Zwg^D. 



-'D{q)6q,q' {ScbSda + ScaSdb) , 

(28) 



(«) 



■■■0 + -0-0 + 



(6) 



FIG. 4. Some Feynman diagrams appearing in the \/Nf 
expansion, (a): summing fermion bubble diagrams [Eq. (Eq)]. 
(b): the 1/Nf expansion correction to the fermion self-energy 
[Eq. (§)]. 



where a,b,c,d~ 1, . 
we have introduced 



,ANr, := {2^f5{q-q'), and 



(29) 



(q2 +4I]2)atan(|g|/2i;)' 

In the UV limit (i.e., at the critical point g = g^), 

I?(g)^4/|q|, (30) 

and hence, by choosing the renormalization constant as 
Zw9^ = 



{Wab{q)Wcd{-q')) 



32 



■5q,q' {^cbSda + Sca^db) ■ 



(31) 



The critical exponent u is given, to leading order in 
1/Nf, by I' = l/{d — 2) where d = 3. We now com- 
pute next to leading order (the order 1/Nf) correction to 
the critical exponent i>. First, we compute the wavefunc- 
tion renormalization and mass renormalization to next- 
to-leading order. The fermion propagator with 1/iV/ cor- 
rections is given by 



g ^{k) = 2Z^YkiliSab 

+ 'SAab 



1 - 



16(4AV + 1) 
iir^NfD^ 1 77 1 

_16 , ^ 

NfD, M 



(32) 



From this, we determine the counter terms as 

A, 



I 



16(4Ar, + 1) ^ 
"f^i \v\ 



E = 77 
with a : 



1 + aln 



A, 



= V 



V 



64(iV^ + 1) 
STT^NfD-y ■ 



(33) 



By definition, the critical exponent u is determined 
from the divergence of the correction length ^ as a func- 
tion of the distance from the critical point. 



^)(xe 



-l/u 



(34) 
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Note also that ^ is inversely proportional to rj. On the 
other hand, to this order, t oc S**^^. Combining these. 



t oc r;(''^2)(l-«) ^ ^-(d-2)(l-a) 

1 _ + 



id -2) {I -a) {d-2)' 



(35) 



In particular, when d = 3, A^/ 
replica limit A'^^ ^ 0, 



1, = 4, and in the 



v = l 



16 
3^ 



1.54. 



(36) 



The result should be compared with the critical expo- 
nent at the conventional Anderson transition in symplec- 
tic (All) symmetry class, which is estimated numerically 
as ^ 1.3 — 1.4. (See Ref. 40 and references therein.) 

a. dynamical exponent So far we have been inter- 
ested in static properties. Some dynamical properties 
can also be computed within the large- A^/ expansion. To 
this end, let us start from the imaginary-time action 



S= drd'^x J2 ^! iidr + no) S,^ + 

t,K = l 



(37) 



+ i 



d'^2 



Within the imaginary-time action, electron-electron in- 
teraction terms such as (A/2) J drd'^x ( V'lV't) can 
be added. As before, disorder can be averaged by 
the fermionic replica method, and then the subsequent, 
disorder-induced "interaction" term can be decoupled in 
terms of the Hubbard-Stratonovich field W{t,t' ^x) as 

S = / drd-'x J2 i9r + no) (38) 

1 

dTdT'J2i^lir)Wabir,T')Mr') 

a,b 

+ Y j '^"^ I dTdT'Y,Wab{r,T\x)WUr',T,x). 

ah 

This action can be studied by the saddle point method. 
We assume the saddle point is spatially homogeneous and 
replica symmetric, 

Wab{T,T',x)=X{T,T')5ab, (39) 

and X{t, t') in frequency space is diagonal 

X{iujn, iLOn') = <5(a;„ - w„')X(iw„). (40) 

The electron Green function is given by, at the 
disorder-dominated critical point, 

Q~^{uin, k) = -iuJn + no{k) + iX{iu!n) 

= Ho(fc)-*C'rf|a;„|^sgna;„, (41) 

where Cd is a dimensionful numerical constant, Cd = 



{[2|sin(7rd/2)|]/[7r(d-2)]} 
this, the dynamical exponent is 

z^d-l. 



l/(d-l) ^(d-2)/(d-l) 



From 



(42) 



The Green's function ( p4| ) should be contrasted with the 
non-rclativistic case where the disorder induced imag- 
inary part of the self-energy within the self-consistent 
Born approximation is given by 



iX{iuJn) = -7;-sgna;„, 
Zt 



(43) 



where r — \/^'np{eF)g\ is the elastic scattering time with 
p(£^ being the density of states at the Fermi energy 



III. DISORDER IN THREE-DIMENSIONAL 
TOPOLOGICAL SUPERCONDUCTORS 

The Bogoliubov-de Gennes (BdG) Hamiltonian de- 
scribing the dynamics of fermionic quasiparticles in the 
superconducting (superfluid) phase realizes, depending 
on its symmetry properties, six out of ten Altland and 
Zirnbauer symmetry classes. They can be summarized, 
in terms of how badly SU(2) spin rotation symmetry is 
broken, as follows: A family (an ensemble) of BdG Hamil- 
tonians with completely broken spin rotation symmetry 
is called symmetry class D (no TRS) or Dili (with TRS). 
When one component (z-component, say) of spin is con- 
served, symipctxv class A (no TRS) and AIII (with TRS) 
are realized. EJO Finally, when full SU(2) symmetry is 
preserved, class C (no TRS) and CI (with TRS) are re- 
alized. 

It is worth mentioning that historically symmetry class 
A and AIII have been discussed mainly in the context 
of electronic systems (non BdG systems): In the origi- 
nal paper by Altland and Zirnbauer, class D, Dili, C, 
and pi are called "BdG classes", but not class A and 
AIII.E3 Symmetry class A, or unitary class, appears quite 
generally for electron systems which lack with TRS, and 
also with any other discrete symmetries. Symmetry class 
AIII describes an electronic system with broken TRS and 
with sublattice symmetry; i.e., tight-binding Hamiltoni- 
ans with bipartite hopping elements only. While sublat- 
tice symmetry implies an invariance of the single-particle 
energy spectrum under reflection about zero energy, this 
is not related to PHS of BdG systems when a class AIII 
system is interpreted as describing BdG quasiparticles. 
Below, when we speak about symmetry class AIII, we will 
focus on its superconducting interpretation, since sublat- 
tice symmetry in electronic systems requires fine-tuning, 
while class AIII in BdG systems arises more naturally. 

In three spatial dimensions, there are three symme- 
try classes, out of ten Altland and Zirnbauer symmetry 
classes, for which non-trivial topological phases charac- 
terized by an integer topological invariant exist. These 
are symmetry class CI, AIII and DHL BdG Hamiltoni- 
ans in these classes have two discrete symmetries, TRS 
and PHS of different kind and as a result they anticom- 
mute with a unitary matrix. For this reason, topological 
phases in these symmetry classes are called chiral topo- 
logical superconductors in Ref. ^6l 
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In this section, we will discuss effects of disorder in 
these 3d topological superconductors. For the most part 
of this section, we will mainly focus on symmetry class 
AIII. The formalism for classes Dili and CI are largely 
parallel to class AIII, and will be discussed in Appendix. 

A. more descriptions of models 

b. BdG Hamiltonian The dynamics of quasiparti- 
cles deep inside a superconducting phase is described by 
a BdG Hamiltonian. On a lattice with N sites, which is 
convenient for the purpose of discussing discrete symme- 
tries of various kind, it can be written as follows: 

H^\{c\ c)'H,^l,y'H,= ^l, _%ry (44) 

where T-L^ is a AN x 4A^ matrix and c — (c-j-, cj,). [c 
and can be either column or row vector depending 
on the context.] Because of S = (hermiticity) and 
A = -A^ (Fermi statistics), the BdG Hamiltonian ( p^ ) 
satisfies PHS 

[PHS] : ^4 = ~t,njt,, (45) 

where the Pauli matrix tx.y.z is acting on the particle- 
hole grading. In addition to PHS, depending on physical 
systems, TRS 

[TRS] : ^4 = iSy-Hji-iSy), (46) 
and rotation symmetry for each spin rotation axis 
[SU(2) symmetry (a = x, y, z)] : 

[n,,Ja] =0, Ja:= _^t), (47) 

can be satisfied, where s^^y^z is the Pauli matrix acting 
on spin indices. Depending on their symmetry properties, 
BdG Hamiltonians (Q) are classified into sub classes D, 
Dili, A, AIII, C, and CI. 

1. class Dili 

When both TRS and PHS are preserved in the BdG 
Hamiltonian, the relevant symmetry class is symmetry 
class DHL Combining TRS and PHS, one can see that 
a member of class Dili anticommutcs with the unitary 
matrix (8) Sy, 

V.^ = —tj. (g) SyT-L/^t^ ® Sy. (48) 

In this sense, class Dili Hamiltonians have a chiral struc- 
ture. It is sometimes convenient to take a basis in which 
the chiral transformation, which is tx ® Sy in the present 
basis, is diagonal. In one of such bases, a class Dili 
Hamiltonian takes on the form 

H4=(^^ o), D = -D^. (49) 



A canonical example of the class DUL-JpnaJogical su- 
perconductor is the B phase of '^Hc.EltjOcJ which is 
described, in momentum space, by the following BdG 
Hamiltonian: 

H=\j d^fc*^(k)H(k)*(k), (50) 

where ^'''(k) = (c| j^., c| C|. Cj^ is the Nambu 
spinor composed of fermionic creation/annihilation op- 
erators {c\ jj./Cj 1^) of a ■^He atom with spin s and mo- 
mentum k, and the kernel takes the following form: 

- ( i.'ri -tf I, ) ^ 

The matrix elements are given by 
e(k) = e/{2m) - ^, 

A(k) = <i(k)-s(is„), d(k) = |A|k, (52) 

where m is the mass of a '^He atom, /i is the chemical 
potential, and |A| is the amplitude of the pair poten- 
tial. With the d- vector pointing parallel to momentum, 
there is an isotropic gap everywhere on the 3d fermi sur- 
face. The critical point at = separates topologically 
trivial {fi < 0) and non-trivial (/x > 0) phases, which 
are characterized by the topological invariant v = Q and 
1/ = 1, respectively. A recent surface transverse acous- 
tic impedance measurement reported in Ref. ^ and ^ 
revealed the existence of the surface Majorana Dirac 
fermion mode on the surface of "^He-B. 

Many recently found non-centiMsymmetric supercon- 
ductors also fall into class DIILE^I The theory of disor- 
dered topological superconductor in class Dili which we 
will present is also relevant to these materials. A Ma- 
jorana hopping model on the diamond lattice, which re- 
alizes trivial {v = 0) and non-trivial (|z/| = 1) phases, 
as well as the Dirac Hamiltonian similar to Eq. ( ]5T| ) was 
discussed in Ref. 

Below, we will consider effects of disorder around this 
transition separating the weak-strong pairing phases, by 
perturbing the BdG Hamiltonian (|^) by inhomogeneous, 
time-reversal invariant potential V{x). With a suitable 
unitary transformation, the full BdG Hamiltonian is writ- 
ten in terms of the gamma matrices (^ as 

n = |A|(-0,)a. + e(-*aO/3 + V{x). (53) 

We will henceforth set |A| = 1 and drop the 0{k^) term 
in C(k), C(k) ^ -Ai. 



2. class AIII 

Let us consider BdG Hamiltonians which are invariant 
under rotations about the z- (or any fixed) axis in spin 
space, [7^4, J^] = 0. This implies that the Hamiltonian 
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can be brought into the form: 

where and S are an NxN matrix and = ^a- With- 
out further constraints, this Hamiltonian is a member of 
class A (unitary symmetry class). If, in addition to con- 
servation of the z-component of spin, we further impose 
TRS, we obtain the conditions = and S = S\ which 
can be summarized as 

ry'H2ry = -7^2, (55) 

where the Pauli matrices r^^y^z are acting on the particle- 
hole grading. It is also convenient to rewrite the Hamil- 
tonian by rotating the r^^y^z matrices by {rx,ry,rz) =4» 
{r^, —rz,ry). In this basis, the class AIII BdG Hamilto- 
nian takes on block off-diagonal form 

^^=(i?to)- (56) 

While we are not aware of any (lattice) BdG Hamil- 
tonian realizing topological superconductor in class AIII, 
an electronic lattice model of class AIII with non-trivial 
topological charge v = 1 is constructed in Ref. In the 
continuum, the model can be described in terms of the 
following massive Dirac Hamiltonian: 

"Ho = -idkctk - ifij^m^. (57) 
3. class CI 

If, in addition to [7^4, Jj,] = 0, we further impose the 
full SU(2) rotation symmetry, and S are constrained 
by = =• ^7 S ^ S'^ . This defines symmetry class C. 
Imposing both full SU(2) rotation and TR symmetries, 
leads, in addition to = and 6 = 5'^ , to £,* = and 
6* = 5. These conditions can be summarized as 

''"y'^I'i'y = ^^2i ^2 = ^2- (58) 

This defines symmetry class CI. Combining these two 
conditions we can obtain a chiral symmetry, ry'H2fy = 
— By rotating the matrices as {rx,ry,rz) ^ 
{r^, —rz,ry), the class CI Hamiltonian takes on block off- 
diagonal form 

^2 = ( i5t 0* ) ' D^S-i^ = D^. (59) 

In three dimensions, class CI is the only class which 
admits non-trivial topological state with SU(2) spin ro- 
tation symmetry. A BCS Hamiltonian on the diamond 
lattice which belongs to symmetry class CI and real- 
izes topological superconductor with i> = 2 was con- 
structed in Ref.pftSf A model on the cubic lattice was 
also constructedtrlTand it was shown that a defect in such 
topological SC accumulate a spin 1/2 degree of freedom 
at its core. 



B. continuum Dirac model and fermionic replica 
method 

We now start our discussion on the effects of disorder 
on topological superconductors. Let us start by consid- 
ering, as a model of a 3d topological superconductor in 
class AIII, a 3d massive Dirac Hamiltonian 

Ho = -idkdk - iPj^ms, 

where Dj (x) in V represents a disorder potential. For any 
value of the mass and for any realization of disorder, 
the Hamiltonian ( ]60| ) satisfies chiral symmetry 

= -H. (61) 

As discussed in Ref. ^ Ho realizes, depending on the 
sign of the mass term 7715 e M, a topologically triv- 
ial and non-trivial superconductor in class AIII, each 
of which is characterized by the vanishing and non- 
vanishing, respectively, of the winding number v de- 
fined for the projector introduced in Ref. |[ We assume 
Di{r) is maximally random according to the distribution 
«cxphtr2(DjL'|)/(2.g)]. 

To discuss effects of disorder, it is convenient to intro- 
duce the generating function, Z ~ J 2?[7/'J, -0^^] exp ( — 
J d^x £) , for the retarded Green's functions, where 
ipg^ with a = 1, • • • ,Nr '\sa, replicated fermionic vari- 
able, and 

C = Y^iPi{i'H + Q+)^l:a. (62) 

a=l 

We have introduced Nr replicas for the quenched disorder 
averaging, with A^^ — > at the end of calculation. The 
O"*" > is introduced for the convergence of the functional 
integral. The Lagrangian enjoys an \][Nr) symmetry 

i/'t^?/'^r^ i^^T^P, T e\J{Nr). (63) 

Due to the chiral symmetry of class AIII, there is an 
additional (axial) \]{Nr) symmetry when O'^ — 0, 

^^^^^e+'^\ i^^e+^'^V, e^n{Nr). (64) 
Averaging over disorder gives 

a 

- I ^ {'4>a1pb'(pb4'a - V'a7°'0fc'/'67°V'a) , (65) 
a, 6 

where = ip^ j3i and 7'"' = —ifia^. [Observe that we 
introduced 7*"' differently from the case of symmetry class 
All, Eq. (p5|). For symmetry class AIII, Dili, and CI, we 
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find this convention is more convenient.] This is a 3d 
version of theuchiral Gross-Neveu or the Nambu-Jona- 
Lacinio modelE3 perturbed by a mass term. 

The quartic interactions among fermions can be decou- 
pled by the Hubbard-Stratonovich transformation with 
two Hcrmitian matrices W and 

a 



a.b 



+ —tvN^ {W^ + V^) 



(66) 



There are similar Dirac type Hamiltonians describing 
3d topological insulators in class CI and Dili, and the 
following mean-field and large- TV/ analysis are completely 
parallel for these classes. In the X&ige-Nf calculations, 
the nature of Nambu-Goldstone fluctuations are different 
among symmetry classes AIII, Dili, and CI: Roughly, one 
needs to replace U — > O for class DIII, and U — Sp for 
class CI. 



C. large A^/ -analysis 

We now look for a spatially homogeneous, replica sym- 
metric saddle point solution, by setting 



Wab = W6ab, Vab = vSab- 



(67) 



It is enough to look for a saddle point with v ^ and w ~ 
0, since there is a symmetry that relates a saddle point 
with w = and v ^ and a saddle point with w ^ and 

V = 0. The resulting mean-field (self-consistent Born) 
phase diagram is essentially the same as Fig. |[ (We 
need to replace m — > ms). We identify the phase with 

V ^ (the ordered phase of the NLctM) as a (diffusive) 
metallic phase. This metallic phase appears because the 
disorder creates states that fill the band gap. While it is 
not depicted here, in the large g limit, there should be 
an Anderson insulating phase in class CI and class DIII. 
On the other hand, in class AIII, it is not necessary to 
have an Anderson insulator in the large g limit since the 
beta function in (2-|-e)|dimensions is positive to all order 
in perturbation theory,E3 as known in the context of the 
random hopping (flux) model in one and two dimensions 
where there is no Anderson localization. For a numerical 
stud y o f the 3d random hopping model in class AIII, See 
Ref. 54, (This result should be considered, however, as 
the case with zero topological angle, — 0, while our 
system is located at 6 = n/2 when ~ 0; see below for 
more discussion.) 

Just at the Dirac point, the disorder, which appears 
as a four fermion interaction in the generating function, 
is irrelevant, whereas the mass term is relevant, from 
the power-counting. The transition between trivial and 
topological insulators is then described by the clean Dirac 
point. 



In the large- iV/ expansion, the model can be renormal- 
ized in a similar way as the class All Gross-Neveu model 
(p5|). We introduce wavefunction renormalizations 
and Zw and consider the Lagrangian 

C = Z^i),a{dklk5ab + ^)i^i.b 

i^ia {Wab + il^Vab) ^Lb 



(68) 




+ l^Zw {WabWab + VabVab) 



where we have rescaled bosonic fields properly. To lead- 
ing order, the fermion propagator and the boson propa- 
gator are given by 

{'4'La{k)ipKbik)) = Si,JabZ:^,^{ikiji + S), 

{Wabiq)W,di-q')) - {Vabiq)Vcd{-q')) 

4 



Zwg'^D, 



-Viq')6g^g,5cb5da. (69) 



As in the class All Gross-Neveu model, several critical 
exponents can be readily computed within the large- iV^ 
approach. For example, the anomalous dimensions of a 
fermion ip to leading order in 1 /Nf is given by 



nip = a 



WNr 



for class AIII, 



(70) 



where the fermion Green's function behaves as Gix) ^ 
|2;|-2-j)v,^ In the replica limit, this leading order correc- 
tion is vanishing. Similarly, the anomalous dimension of 
the fermion mass is given by 



32Nr 



(71) 



and the anomalous dimension of the "75" fermion mass 
is given by 



32Nr 



(72) 



The similar calculations go through for class DIII, and 
CI. For symmetry class DIII, 



64iV, 



(73) 



and for the anomalous dimension of the fermion mass is 
given by 



256 {Nr - 3) 



(74) 



D. non-linear sigma models 

Deep in the diffusive metallic phase, physics at large 
distances should be described by the NLcrM defined on 
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symmetry 


NLcrM 


3d topological 


class 


target space 


term 


A 


U(2iV)/U(iV) X U(iV) 




AI 


Sp(2iV)/Sp(7V) X Sp(iV) 




All 


0{2N)/0{N) X 0(iV) 




AIII 


U(iV) X U(iV)/U(iV) 


Z 


BDI 


U(2iV)/Sp(iV) 




CII 


U(27V)/0(2iV) 




D 


0{2N)/U{N) 




C 


Sp(iV)/U(iV) 




Dili 


0{N) X 0(iV)/0(iV) 




CI 


Sp(iV) X Sp(Af)/Sp(iV) 





TABLE I. Fermionic replica NLcrM target spaces and possible 
topological terms in three dimensions. For a similar table for 
two spatial dimensions, see Ref. M . 



U(iV^) (AIII), 0{2Nr) (Dili), and Sp(7V^) (CI), respec- 
tively. The coupling constant of the NLcrM plays the role 
of the longitudinal conductivity Oxx- (To be more pre- 
cise, in superconductors, Oxx corresponds to either spin 
or thermal conductivity.) In addition, one should note 
here that the NLcrM for these classes AIII, Dili, and CI 
can support a topological term in three dimensions: For 
the NLcrM target manifolds G, where G = U(iVr) (AIII), 
G = Q){2Nr) (Dili), and G = Sp(iVr) (CI), the homotopy 
group is given by 7r3(G) = Z. 

The NLcrM and the topological term can be obtained 
in the following way (we will focus on symmetry class 
Am as an example, but a similar derivation should be 
possible for class Dili and CI): a saddle point solution 
in Eq. ( |67| ) breaks the global U(Afr) x U(iVr) symmetry 
down to ViiNr). One then expects the long wave length 
physics are described in terms of the Nambu-Goldstone 
modes associated with the symmetry breaking. With the 
inclusion of the Nambu-Goldstone modes (=fluctuations 
around the saddle point), the action is 



where the slow variations represented by exp(2^7*'^^(x)) 
are the Nambu-Goldstone modes in the NLcrM. 



U 

W 



U e \J{Nr). 



(76) 



The Lagrangian ([75|) is obtained from Eq. (66) by freez- 
ing the longitudinal fluctuations but keeping the trans- 
verse fluctuations (i.e., the Nambu-Goldstone modes). 
The longitudinal fluctuations are gapped and can safely 
be integrated over, which can renormalize parameters in 
the effective action. (In symmetryj-dass AIII, it can also 
generate the so-called Gade term.E3 We do not discuss, 
however, the effect of the Gade term in this section.) The 
effective action S'cff[?7] for the Nambu-Goldstone modes 
is then obtained by integrating over fermions. 



V [U] e 



-Sett [(7] 



(77) 



Expanding the effective action in powers of U and its 
gradient, and with the Pauli-Villars regularization£3 the 
effective action (the NLcrM action) is given by 



(fxti {djU^'djU) + ier[U], (78) 



where the topological term 



nU] = j ^e'ihY{U^d,UU^d,UU^dkU) (79) 

is an integer for any field configuration. The 6 angle 
(topological angle) is given by (see Appendix ^) 



\M\- 



\mR\ 



4|M|3 ' 



:= ml 



(80) 



where mg is the imaginary part of the self-energy (within 
the self-consistent Born approximation), and is a 
mass introduced by the Pauli-Villars regularization. No- 
tice that the topological angle is tt when is (Dirac 
point) for arbitrary mo. 



1. possible two-parameter scaling phase diagram 

The presence of the topological term is a direct mani- 
festation of the fact that the metallic phase is adjacent to 
a topological superconducting phase. The similar topo- 
logical term in two dimensions, the Pruisken term in the 
QHE, has a significant implication on the nature of An- 
derson localization physics, and, in fact, is crucial for the 
existence of the QHE, and for the plateau transition. The 
role of such Pruisken term in topological superconductors 
with broken time-reversal symmetry in two dimenskma 
has been discussed in symmetry class A, D, and C.EZrEa 

It is not clear what the role is played by the topolog- 
ical term in 3d NLcrMs. While controlled calculations 
near metal-insulator transitions and in the presence of 
topological terms are in general difficult, below we com- 
bine the RG flow of the NLcrM for large axx with the 
RG flow around the clean Dirac point from the l/Nf- 
expansiO|ii, to speculate the general structure of the phase 
diagramed As in the QHE, the phase diagram and the 
RG flow may be phrased in terms of the two parameters 
of the 3d NLcrM; the NLerM coupling constant 1/(Txx and 
the topological angle 9x3 

We first discuss the phase diagram for class CI and 
DHL (Symmetry class AIII has a complication due to the 
special nature of its beta function in (2 -I- e) dimensions 
- see below). The beta functions for these NLcrMs in 
3d when 9 = are expected to have a zero at finite 
axx, separating stable metallic and Anderson insulating 
phases. For large axx, the RG flow for arbitrary 9 is 
expected to be similar to the RG flow at 6* = 0, while the 
flow for 9 and for small axx needs to be determined. 
In terms of the massive Dirac model with disorder, on the 
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277 TT 2tT 



FIG. 5. Possible two-parameter scaling flow diagrams of 
symmetry class CI and DHL 



other hand, axx is smah around the clean Dirac fermion 
point (the dc conductivity of the clean 3d Dirac fermion 
is zero) and 9 ^ tt. It is thus tempting to "identify" the 
point {(TxxtO) = (0, tt) in the NLcrM phase diagram as 
the clean Dirac fixed point, or, at least, to use the Dirac 
model (i.e., the chiral Gross- Neveu model) to guess the 
RG flow of the NLcrM. 

The fermion mass term ip^^ip in the chiral Gross-Nevcu 
model is a relevant perturbation to the clean Dirac fixed 
point. Also, depending on the sign of the mass, one can 
drive the clean Dirac fermion either into the topolog- 
ical {f ^ 0) or trivial {i' = 0) insulating phase. As 
the fermion mass corresponds to the theta term in the 
NLcrM, this would suggest the deviation of the theta an- 
gle from 9 = TT is relevant at a^x = 0. On the other 
hand, disorder at the clean Dirac fixed point is irrele- 
vant by power-counting. Furthermore, the existence of 
the UV fixed point (the Gross-Neveu fixed point) in the 
l/A^^j-expansion of the chiral Gross-Neveu model would 
suggest the existence of a corresponding critical point at 
9 = n ill the NLcrM (Fig. ||) separating a metallic phase 
and an insulating phase (which is adiabatically connected 
the clean 3d Dirac fermion with axx ~ 0). 

While we expect the deviation of the theta angle from 
= TT is relevant at cFxx = (i.e., at the clean Dirac fixed 
point), there is an ambiguity on the RG flow around the 
putative critical point at = tt and axx 7^ 0. If it is 
relevant, as it is at the clean Dirac point {axx, 9) = (0, tt), 
the flow looks like Fig. ||-(a). An alternative scenario 
where the theta term is irrelevant at the Gross-Neveu 
critical point is depicted in Fig. ||-(b). In these scenarios, 
an implicit assumption we have made is that there are 
no other fixed points. Also, we have not considered more 
exotic possibilities, such as the case where 9 is exactly 
marginal at the Gross-Nevcu fixed point. Which one of 
the two scenarios is realized may be speculated, again, 
from the large Nf approach, by computing the anomalous 
dimension of the 75 fermion mass; See Eqs. (^, (|72]), 
and (^Sj) for relevant calculation. We should, of course, 
bear in mind that the structure of the phase diagram 
would be quite different for large Nf and for small Nf 
(e.g., Nf — >■ 1). In other words, with the inclusion of 



the Nf flavor of Dirac fermions, the model describes, in 
the absence of disorder, a transition between trivial and 
topological insulators with v = Nf. This is a highly fine 
tuned critical point, and might be quite different from 
the critical point that may exist and be described in the 
NLcrM. 

Symmetry class AIII is special in that the beta func- 
tion in (2 -|- e)-expansion is always positive; Thus, at 
9 = 0, there is no metal insulator transition, but we 
always have a metal. However, as the large-iV^ expan- 
sion (which corresponds to 7^ 0) predicts tuning 6* 7^ 
seems to generate a metal-insulator transition. Because 
of the positivity of the beta function in (2+e) dimensions, 
it appears that constructing a reasonable conjecture on 
the two-parameter scaling phase diagram is more com- 
plicated. 



IV. DISCUSSION 

Topological insulators and superconductors are known 
to be an inherently holographic system; The non-trivial 
bulk topology [as encoded by a non-vanishiim-bulk topo- 
logical invariant such as the Chern (TKNN)Ej integer or 
the Z2 invariant in Z2 topological insulators] has a crucial 
impact on the Anderson localization problem at bound- 
aries of the system; in fact, Anderson localization is not 
allowed at the system's boundaries because of topological 
reasons.EaEa Such bulk-boundary correspondence served 
as a guiding principle of classification of topological in- 
sulators and superconductors. In this paper, we have 
investigated Anderson localization problem in the bulk of 
topological insulators and superconductors in three di- 
mensions; we have explored implications of bulk topology 
on physics of Anderson localization in the bulk. 

In the QHE, the transition between topological phases 
with different Chern integers is governed by a disorder- 
dominated quantum critical point (the plateau transi- 
tion). The fact that delocalized (extended) states are 
possible at the plateau transition is a direct manifesta- 
tion of a topological character of the system, which is not 
allowed otherwise in two dimensions with broken time- 
reversal symmetry. In the NLcrM description of the QHE, 
the plateau transition is encoded in the Pruisken term 
(which also points the existence of edge modes). We 
can draw a rather close analogy between time-reversal 
symmetric topological superconductors in three dimen- 
sions and the 2d QHE, both of which are characterized 
by an integer topological invariant in the bulk. We have 
given a field theory description of metal-insulator transi- 
tion in terms of the Dirac model of topological supercon- 
ductors by employing the large- A^j expansion, and also 
in terms of the NLctMs. Although the presence of topo- 
logical terms in the NLcrMs is indicative, it is an open 
question if the phase diagram is characterized in terms 
of two parameter scaling, as in the QHE. It is desirable 
to explore further the structure of the phase diagram and 
critical behaviors in 3d topological superconductor sys- 
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terns in terms of, say, numerics, preferably in the presence 
of the non-zero theta angle in the language of the NLcrM. 
This would require more microscopic and physical under- 
standing of the theta angle; e.g., how one should iptetpret 
and measure the theta term in a metallic phase.r^rl 

Our discussion for symmetry class Dili can be relevant, 
for example, for "^He in aerogel,E§E3 although we need to 
include pair breaking effects. And-jalsp Cu-intercalated 
topological insulators, Cu2:Bi2Se3,E2rO and many heavy 
fermion non-centrosymmetric supfirconductors with an- 
tisymmetric spin-orbit coupling£3 Experimental probes 
for detecting these systems are, for example, thermal 
transport, and spin transport phep spin or part of spin 
quantum number is conserved. EilO See, for example, for 
the anomalous thermal Hall effect by magnons, Ref. |7^. 

In the 3d Z2 topological insulator in symplectic sym- 
metry class, a closer analogy would be the 2d quan- 
tum spin Hall effect. We have numerically investigated 
the structure of the phase diagram by the transfer ma- 
trix method, and developed the large-iV/ field theory 
approach to a metal-insulator transition. It should be 
noted, however, the large Nf expansion does not distin- 
guish the case of odd Nf , for which the one of the insulat- 
ing states in the phase diagram Fig. ^ is a topologically 
non-trivial, Z2 topological insulator, and the case of even 
Nf, for which the one of the insulating states in the phase 
diagram Fig. || is a topologically trivial. While we have 
imposed the global 0{Nf) symmetry in the large Nf ex- 
pansion, perhaps a distinction between even and odd Nf 
would emerge once we study the stability of the UV fixed 
point (the Gross-Neveu fixed point) that can be reached 
by the large Nf expansion; for example, for even Nf, 
the Gross-Neveu fixed point is unstable against 0{Nf) 
symmetry breaking (and other) perturbations, and the 
metal-insulator transition is always controlled by the con- 
ventional critical point of symplectic symmetry class. In 
terms of the NLcrM description, there is no topological 
term of any kind in the bulk in symmetry class AH, and 
hence the transition between the metal and insulating 
phases is expected in the conventional universality class 
of symplectic symmetry class, irrespective of the number 
of flavors. This situation is similar to thc_Zg^ topological 
insulators in class AH in two dimensionsOu. 

While we have discussed the bulk critical behavior, it 
would be interesting to study boundary critical behav- 
iors; even when a transition between a trivial insulator 
to a metal and one between a topological insulator to a 
metal is in the same universality class, at a surface we 
might observe different critical behavior. This was indeed 
the case in the 2d symplectic class.LZl 

Finally, we conclude with some remarks on other sym- 
metry classes in three dimensions. 

- For symmetry class CIL-|a Z2 topological insulator 
can be realized in the bulk.Ej While we did not discuss 
symmetry class CII, a similar field theory analysis can be 
developed. 

- The large-TV/ technique developed here is also ap- 
plicable to the Weyl semimetal in symmetry class A, 



which arc .realized in the A phase of ^He and Pyrochlore 
Iridates.EaOl 

- We mention two other symmetry classes for which 
the corresponding NLcrM has a topological term of Z2 
type: in symmetry class C and CII, Z2 insulators (super- 
conductors) can be realized in four dimensions. Thus, 
the 3d Z2 topological term in the corresponding NLcrM 
describes the surface of 4d insulator, but not any isolated 
3d system. 
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Appendix A: transfer matrix 

a. transfer matrix In this appendix, we derive the 
transfer matrix for four-component 3d Dirac fermions. 
The Dirac Hamiltonian is written in the single-particle 
representation as 

^ = - II [1^' " + 1' ^)"-,aa' (r, n, a'\ 

n.r 

+ \r,n,a)ul,^^,{r,n + l,a\] 

[\T^ + et„n,a)u'^^^^,{r,n,a'\ 

n.r fi— 1.2 

+ |r,7i,cr)u^,^^,(r + e^,ri.,cr' 
+ |r, n, cr) [(to + 3)/3 + ern]aa' (r, n, cr'l , (Al) 



where we have introduced 



(A2) 



The eigen wave function with energy E of this Hamilto- 
nian satisfies 

0== {rna\{H - E)\'iP) 
= -ut^„„4'a' (r, n - 1) - u% ,^„,i}}„, (r, n + 1) 
+ [H(")-i?],,-^v(r,n), (A3) 

which leads to 

V'n+i = 2(/3 + ^^a,)-M^'"' - -ElV'n 

~ + ia,)~\p - ia,)i^^-i. (A4) 
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Note that (P + ita,)-'^ = {l3 + ita^)/{l-t'^), and Eq. 
is obtained. The transfer matrix M'^"^ is defined as 



where 



= X! |r + e^,n,cr) 

n,r,/^— 1,2 



4ui?il 



t2- 1 



n,r,/^— 1,2 



4lt^M^ 



(r, n, (7 



t2- 1 

^|r,n,t7)^^[(m + 3)u;/3 

+ (er„-i;Kl {r,n,<j'\. (A6) 



Appendix B: large- A*'/ expansion for 
three-dimensional topological superconductors 

In this Appendix, we ihustrate the large- A^/ calcula- 
tions for 3d time-reversal symmetric topological super- 
conductors in s ymm etry class AIII, Dili, and CI. We 
start, in Subscc. B 1 , by reviewing fermionic replica meth- 



ods for sy mmet ry class AIII, Dili, and CI. The formalism 
in Subscc. B 1, is general and applies to any systems in 
these symmetry classes in any dimensions other than 3d 
Dirac representatives of topological su perc onductors dis- 



cussed in the main text. In Subscc. 



B2 



we introduce 



auxiliary matrix fields (Hubbard-Stratonovich fields) to 
decouple four fermion interactions induced by quenched 
disorder averaging, which is useful to develop the large- 
Nf expansion, and also to derive the NLcrM inside metal- 
lic phases. 



1. fermionic replica field theories 

For symmetry classes AIII, Dili, and CI, the Hamilto- 
nians anticommute with a unitary matrix with equal 
number of eigenvalues -fl and —1, {HjTz} = 0, thereby 
they can be brought into an block off-diagonal form, 

^-(7?t?)> (Bl) 
in the basis where is written diagonally as 

For class AIII, the Hamiltonians arc not subjected to any 
further constraint. On the other hand, for classes Dili 
and CI, the off-diagonal block D satisfies 



= -D for class Dili, 
= +D for class CI. 



(B3) 



In the following, we will assume that the Hamilto- 
nian consists of the disorder-free (Ho) and disordered (V) 
parts. 



n = no + V, D^Do+Di{x), 



(B4) 



where the random potential V is assumed to be white- 
noise and to be maximally random according to the dis- 
tribution 



(X exp 



-ti-{DjDl)/{2g) 



(B5) 



with g characterizing the strength of disorder. 

In order to study Anderson localization physics, we are 
after the properties of the single particle Green's func- 
tion G{z) = {z — H)^^ at zero energy, Rez = 0. Since 
"chiral" symmetry {7i,r.} ~ relates retarded and ad- 
vanced Green's functions at zero energy, we can focus 
either one of them, the retarded Green's function, say, 
G{z = iO+) = (iO+ - n)-\ where 0+ > is positive 
infinitesimal. The generating function for the retarded 
Green function can be expressed in terms of fermionic 
functional integral, Z = J /] cxp(— / d^x C), 



c^p (in + 0+) /, 



(B6) 



where and / arc independent fermionic path integral 
variables. 

To perform the 1/iV^ expansion, we generalize the 
above functional integral with single flavor to Z = 
/P[/t,/Jexp(-/d3x£), 

= E [(*^" + + *"^-] (B7) 

l,K—l 



where the disorder part 







(B8) 



is a suitable generalization of Hi which couples all flavors 
equally according to the disorder distribution oc exp [ — 

tr(^/,..i?Lj/(25)]. 



a. class AIII 



Quenched averaging over disorder Dj can be done by 
introducing replicas {/^^, fai}a=i,...,Nr which leads to the 
replicated Lagrangian 



— — ( i/t*!" t/i- t/i"^ 



(B9) 



where we have renamed f\f — ipK'^'i and repeated 
replica and flavor indices are implicitly summed. 
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When we deal with the relativistic model (massive 
Dirac model) of the topological superconductor, it will 
prove convenient to introduce 

?A:=V^«r^- (BIO) 
The Lagrangian is then given as 

^Ani=V5a,,(r.Ho-ir,0+)Va. (Bll) 

This is the chiral Gross-Neveu model with U(-/Vr) xU(A^r) 
internal symmetry. 

b. class Din 

For class Dili and (and class CI, which will be dis- 
cussed later), it is advantageous to double the number of 
components of fermionic fields to account for the sym- 
metry of the Hamiltonian, D = ±D^ . For class Dili, we 
first note that 

J d'xfiDfs = + J d'xflDn, 

I d'xflD^f^ = +l d^xf^D^Ts, (B12) 

where / = (/^, Z^)^ hi the chiral (or "sublattice" ) grad- 
ing. Introducing the "charge-conjugation" space by 

the Lagrangian can be written as 

-Cdiii = X+'i'Tx ® Dx+ + X^iTx <E) D^x^ 
0+ 

+ -J {xItxT.X- + X^r,T,x+) , (B14) 

where the Pauli matrices To^x,y,z acts on the "charge- 
conjugation" space. Further redefining the fermionic 
fields as 

71+:=T^yX+, V~-=TzyX-, (B15) 

where r^y = (r^ + Ty) /V2, rJ^T^y ^ -ir^, 
the Lagrangian can be written as 

0+ 

+ -Y + , (B16) 

where we noted TzyT^Tzy = Ty. 

At this stage, we notice that f^Hf is invariant under 
0(2) X 0(2) transformations, 

r?+->0+??+, V-^0_7j_, (B17) 



where 0± G 0(2). On the other hand, the "smearing 
term" /^O"*"/ breaks it down to 0(2). Indeed, in order 
for 0+ and 0_ to leave the smearing term invariant, 
they cannot be independent and must be "locked" by 
0_ = —iTyO^iTy, where iTy S 0(2). With replicas, the 
0(2) X 0(2) symmetry is promoted to 0{2Nr) x 0{2Nr) 
symmetry which is broken down to 0(2Nr) by the smear- 
ing (see below). 

Introducing the notation, 

^t=^(^T ^T)^ (B18) 

the Lagrangian is written as 

Cmu = ip'' {iTo ®na + 0+Ty ® /) V, (B19) 

where we have rescaled 0+/2 — > 0+. Observe that tp^ and 
■0 satisfy the "Major ana" condition 

ij^iV^ = (B20) 

With Nf flavors and Nr replicas, the Lagrangian after 
quenched disorder averaging can be cast into the form 
similar to the one for class AIII, 

-Cdiii - i^l \ino6ab + 0+{Ty (E) /wjah] A, (B21) 

where a, b represents the combined charge conjugation 
and replica index running from 1 to 2A^r- 
If we introduce 

iP := V^»r^, V'^zF.y = (B22) 

the Lagrangian is then given as 

£0111 = ^a, [^znoSab + 0+T,i-lTy)at] A, (B23) 

This is the chiral Gross-Neveu model with 0(2iVr) x 
0(27VV) internal symmetry. 

c. class CI 

To account for the class CI symmetry = +D, we 
note that 

J d^xf\Dfg = - j d^xflDfl, 

I d^xflDUA = - j d^xflD^TB. (B24) 
and introduce 
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The Lagrangian is then written as 

Cci - xlry (E) Dx+ + X-Ty <E) D^^X- 

+ Y [xlHry)x. + X-Hry)x+] ■ (B26) 

At this stage, one notices that f^Hf is invariant under 
transformations 



x+ X- o_x- 

where satisfy 

O^lTyO^ = ITy. 



(B27) 



(B28) 



I.e., 0_|_ S Sp(l). The system is invariant under continu- 
ous Sp(l) X Sp(l) rotations. The smearing term breaks 
this symmetry down to its diagonal subgroup; in order 
for 0± to leave the smearing term invariant. 



O^lTyO- = ITy => 



0+ = 



With replicas, the Sp(l) x Sp(l) symmetry is promoted 
to Sp(AV) X Sp(AV) symmetry which is broken down to 
Sp(AV) by the smearing term. 
Let us now introduce 



^^ = (x?Hr,), X^Hr,)), 



(B30) 



They satisfy the "Majorana" condition 

= V-^r.Hrj,). (B31) 
The Lagrangian is written as 

Co, = V-^ (to ® ino + 0+)^=Y,^l (iHo + 0+) Vr, 

(B32) 



where in the second line we made the particle-hole indices 
explicit, and we have rescaled 0+/2 — 0+. 

Introducing Nf flavors and Nr replicas and then per- 
forming the quench disorder averaging. 



^Cl = V'L,, (^■^O + 0+) S^rSabArL 



(B33) 



where a, 6 = 1, ... , N^, l, k = 1, . . . , Nf, and cr, r = ±. 
If we introduce 

^b:=^phT„ ^p'^{iTy)ry=i;, (B34) 

the Lagrangian is then given as 

^c: - (r.Ho - 0+^^,) (B35) 

This is the chiral Gross- Neveu model with Sp(A^r) x 
Sp(iVr) internal symmetry. 



2. auxiliary matrix fields 

The four-fcrmion interactions in the replica field theo- 
ries can be decoupled by two auxiliary matrix fields which 
we call W and V . Depending on symmetry class, these 
matrix fields satisiy different constraints. 



a. class AIII 

For symmetry class AIII, the fermionic Lagrangian 
(Bll) can be rewritten 



+ i'a.L {Wab + i^zVab) tpb,. 



(B36) 



(B29) where W and V are an Nr x N^ hermitian matrix. 



6. class Din 

For symmetry class Dili, the fermionic Lagrangian 
(B23) can be rewritten 



+ V^a. {Wab + iTzVab)Ac 



(B37) 



where is a 2Nr x 2Nr real symmetric matrix, while V 
is a 2Nr X 2Nr pure imaginary antisymmetric matrix. 



c. class CI 

To decoup le th e four fcrmion terms in the fermionic 
Lagrangian ( Bll ) for symmetry class CI, we need to 
introduce two bosonic fields War.br' and Var.br' where 
the first type of index a, b runs over the replica grading, 
a, 6 = 1, • • • , iVr, and the second type of index r, t' runs 
over the particle- hole grading, r, r' = ±. Wi th these 
auxiliary fields, the fermionic Lagrangian (Bll) can be 
written as 



+ i'aai {Waa,bT + i^zVaa^br) An 

+ ^tr2 (n^Kb + W,,W,,) . 



(B38) 



Here in the last line, the 2d trace is taken over the 
particle- hole space, and Vab and Wab, for the fixed replica 
indices a and &, are viewed as a 2 x 2 matrix. They 
are a real quartenion matrix, i.e., they can be expanded 

as Vab = A + B^{iTx) + By{iTy) + Bz{iT^) and Wab = 
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C + Dx{iTx) + Dy{iTy) + Dz{iTz) with real coefBcients 
Bx.y.zi C, D^yz, and also satisfy 

Ha = (B39) 

3. leading order analysis 

To perform the 1/-/V/ expansion, we first compute the 
propagators of W and V bosons dressed by fermion bub- 
bles. This is depicted diagrammatically in Fig. ^. 

For symmetry class AIII, propagators for W and V 
bosons are given by 

{Wab{q)Wcd{-q')) = {Vab{q)V,d{-q')) 
4 



NfD 



V{q)6g,g<6cb6da, (B40) 



where is defined in Eg. (^. Similarly, for class 

Dili, 

{Wab{-q)Wcd{q')) = ■T^'Diq) [ScbSda + Sca5db) Sq.q, , 
{Vab{-q)Vcd{q')) = ,j „ ^(g) ('^cb(5da - ScaSdb) Sq.q' , 

(B41) 

for a,b,c,d — 1, • . . , 2Nr, and for symmetry class CI, 



-T^iq)K,' 



{WaaM~q)Wca',dr'{q')) - ^^^^ 

{VaaM-q)Vc.',dr'{q')) = 



(B42) 



U(7Vr)yi symmetry. Such continuous rotation is nothing 
but the Nambu-Goldstone mode which is the diffuson or 
Cooperon representing the diffusive motion of fermionic 
quasiparticles. The coupling of the Nambu-Goldstone 
mode to fermions are described by 

£ = {i;{'y''dk + mj"^ + vj'^e^''^"'^^''^)il; + const., (C2) 

where, with space-dependent 0, e'xjp(2ij^9(r)) repre- 
sents the Nambu-Goldstone modes, i.e., the fields in the 
NLcrM, 



^ " ' c/t 



U e V{Nr). 



(C3) 



In the following, we will integrate over fermions to de- 
rive the effective action of the Nambu-Goldstone mode. 
We will do not include the longitudinal fiuctuations for 
a while, but it is impor t a j at , as the integration over them 
leads to the Gade term.E30. To prepare for the gradient 
expansion, we will "unwind" U : 



£ = i: 



mo 



o-kdk + m5 



U {-<Jkdk + TO5) 



mo 



-mo 

o-kdk + mz 

-mo 



^ (C4) 



where we have renamed v — > mo and m — > 771,5. Thus, 
our action can be written as 

£ = ^ {j^dk + mo7° -f 771575 4- -i^Ak -f l^l^Bk) 
where Ak = -Bk = {l/2)UdkUK (C5) 

The effective action is obtained by integrating over 
fermions, 



-Jd^'xC 



(C6) 



Appendix C: non-linear sigma models and 
topological angle for three-dimensional topological 
superconductors 

In this appendix, we will derive the non-linear sigma 
model with topological term for symmetry class AIII. A 
similar derivation should be possible for symmetry class 
Dili and CI as well. 

We start from the saddle point solution Ec^. (|67|): 



{w,v) = (0,w) 



(CI) 



where we have set w = since such saddle points can be 
reached from (0,?;) by continuous \J{Nr) x \J{Nr)A rota- 
tion. I.e., from this saddle point, one can generate other 
saddle points as 717" exp(2i7°6') where 9 G u(7Vr). In 
other words, non-zero saddle point v destroys the "axial" 



and it can subsequently be expanded in powers of Ai, Bi 
and their gradients as 

S-eff = -Tr (Go 1 - V) 

= -Tr In G^^ + V -Tr [GoVf , (C7) 



where V = 'f^Ai + 7^7' S,;, and the free propagator is 
given by 

^ f,\ j'^ikk + mo-f° + m5-f^ 
Go(K) „ , ^2 : ^2 ' i'^^) 



+ TTig -I- m| 



with = mp -I- to| . 



To the 3rd order in the gradient expansion, the effective 
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action is [S^f = (l/n)Tr (Gq^)"] 



Auv 

3tt 



^tr 



'off 



67r|M|' 



'5 Ajk 



47r|A/|3 
By noting that 



e'^'^tr [A.AjBk) 



(C9) 



y d^^e'^'^trA.ajBfc = -2 y dPxe'^^ii B.BjBk, (CIO) 
the combined 2nd and 3rd order terms can be written as 

- [^m,ml+Amf\ 2^iT[U]. (Cll) 

Observe that when mo — 0, if we introduce new gauge 
fields A± by 

A^\{A++A-), B = ]^{A+-A-), (C12) 
the effective action can be written as the double Chern- 



Simons theory: 



5eff = +^ j dPxiv {AjA, + BjB, 



(C13) 



1 ms i 

2 ItosI 47r 



E 

(T = ± 



X tr ( A^d.Al + =-A^A^A^ 



a J d^xe'^'' 

2 

3' 



The effective action ( pll ) [and ( pl4 ) when mo = 0] 
that we have derived should be regularized as signaled by 
the appearance of the term linear in Auv in (Cll) and 
(|C14| ). We introduce the regularized action S^^^y the 
Pauli-Villars regularization with the regulator mass mji 
as 

S'^g-' = 5'off(m5,mo) - lim S'off (mj?, 0). (C14) 

|m_R|-i.oo 

We choose the sign of the regulator mass as 

sgn(mi?) = -sgn(m5) (C15) 

in which case ms > is a topological/non-topological 
phase. We then conclude 



'eff 



+ 



87r|M| 
-1 , 



IMI 



Omsmg + 4m^) 



2\mR{ 



(C16) 
27nr[U]. 
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